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Abstract. We review the status of theoretical evaluations of heavy quark and heavy 
quarkonium hadroproduction cross sections and their comparisons with experimental 
data. 



INTRODUCTION 

When, more than twenty years ago, charm was discovered [1] and subsequently 
interpreted [2] as the first heavy quark ever observed, it came as a big surprise. 
The discovery of bottom [3] a few years later produced perhaps (I wasn't around 
at that time...) less excitement, but was surely equally important in extending our 
knowledge of what is now called the heavy quark sector. 

This sector is today well known, and probably recently completed by the discov- 
ery of the top quark [4] . Theoretical physics of heavy quarks has therefore shifted 
gear, and moved from "discovery mode" to "precision physics." 

The new name of the game is now testing Quantum Chromodynamics, by check- 
ing its predictions against experimental results. The latter have now grown quite 
accurate, and therefore demand equally precise theoretical calculations. In this talk 
I shall describe the state of the art of such calculations for heavy quark hadropro- 
duction. 

I shall first review the fixed-order next-to-leading order (NLO) QCD calculation, 
now available for total cross sections and one- and two-particles distributions. It is 
a consolidated result and provides a benchmark for future developments. 

Large logarithms appear in this calculation and potentially make it less reliable 
in some regimes: log(S'/m 2 ) and log(p^/m 2 ) become large when the center of mass 
energy \/S or the transverse momentum px of the observed quark is much larger 
than its mass. Large log(l — Am 2 /s) appear when the heavy quarks are produced 
close to the partonic threshold. I shall briefly describe the resummation of log(l — 
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4m 2 /s) and \og{p^/m 2 ) terms, and the inclusion of non-perturbative fragmentation 
effects. 

I shall finally also briefly comment on the subject of heavy quarkonium produc- 
tion, where our understanding seems to have been greatly increased by a lot of 
recent theoretical activity. 



I NLO CALCULATION 



The road to the NLO evaluation of heavy quark hadroproduction cross sections 
was paved by Collins, Soper and Sterman [5], who argued that the following fac- 
torization formula holds: 

da{H x H 2 -^QQ ;m ) = J2 I fi/Hj j/H2 da{ij - QQ; m) + O (^9£R\ , (i) 

ij \ m / 

where the summation over partons i and j runs only over gluons and light quarks, 
and the heavy quarks are generated only at the perturbative level, by gluon split- 
ting. The cross section explicitly depends on the heavy quark mass m and on 
all other scales entering the problem (total energy \/S, transverse momentum p?, 
etc.). 

Along these lines, explicit calculations were performed by two groups, Nason, 
Dawson and Ellis [6] on one side and Beenakker, Kuijf, Meng, van Neerven, Schuler 
and Smith [7] on the other. More recently Mangano, Nason, and Ridolfi (MNR) [8] 
have presented a Montecarlo integrator, based on the first of the two calculations, 
which can provide fully exclusive cross sections, thereby allowing detailed compar- 
isons with experimental data. A very extensive collection of such comparisons is 
presented in a recent review [9], from which we select some plots to be shown here. 

The ratios of next-to-leading over leading order predictions for total cross sections 
depend on m/y/S and are about 1.3 for top production at Tevatron energy (yfS = 
1800 GeV) and of order 2 or larger for charm and bottom already at fixed target 
energy. Large uncertainties, due to monotonic renormalization/factorization scale 
dependence, are present for charm and bottom, while the prediction is fairly reliable 
for top (±10%), as shown in figure 1. 

One-particle inclusive differential cross sections, usually px distributions, can also 
be considered. Fig. 2 shows on the left a comparison of E769 pion-nucleon fixed 
target data with NLO QCD plus two non-perturbative contributions: an intrinsic 
initial transverse momentum kr of the colliding partons, with (kr) = 1 GeV, and 
fragmentation effects of the charm into the observed charmed hadrons, described 
with the aid of a Peterson [10] fragmentation function (FF) with e = 0.06. The 
overall result (the dot-dashed line) can be seen to be slightly softer than the data. 
Uncertainties are however large: a larger mass, a larger (kr) or a harder FF could 
help improve the agreement. We should mention that, while the inclusion of the 
non-perturbative contributions might appear unnecessary here, it looks however 
mandatory when considering two-particle distributions or even p T distributions in 
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FIGURE 1. Charm and bottom production at fixed target experiments vs. NLO theoretical 
predictions (left), and top production at the Tevatron (right), from rcf. [9]. 
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FIGURE 2. Charm pt inclusive distribution in nN collisions vs. theoretical predictions [NLO 
calculation + initial kr ((k^) — 1 GeV 2 ) + non-pert, fragmentation (Peterson, e = 0.06)] on the 
left, and azimuthal correlation of the cc pair on the right. From ref. [9]. 



photon-hadron collisions, pure NLO QCD being there clearly unable to describe the 
data by itself. An example is given in fig. 2, on the right, where WA92 pion-nucleon 
data for the azimuthal distance A<p of the c and the c are compared to theoretical 
predictions: the inclusion of the kr kick brings the curve into agreement with the 
data. More such comparisons can be found in [9]. The overall picture suggests 
some consistency between the non-perturbative inputs needed to describe the one- 
and the two-particle distributions in both photo- and hadroproduction. 

More comparisons of pt distributions can be done with the data for bottom pro- 
duction taken at the Tevatron by the CDF and DO experiments. These distributions 
initially caused some concern, as they were markedly higher than the NLO QCD 
predictions. The data have now come down a bit, thanks to a better understanding 
of some decay chains used in the analysis (like the B — > J/ip — > one) and 

to the use of microvertex silicon detectors which allow a much better identification 
and reconstruction of the heavy quark events. But, still, theory and data are not 




FIGURE 3. Bottom pr distributions data versus theoretical predictions, from ref. [9]. 



in perfect agreement. Fig. 3 shows on the left the data from CDF and on the right 
the data/theory ratios from CDF, DO and UAl. Data can be seen to overshoot the 
central theoretical predictions by factors between two and three, and to be in fair 
agreement with the upper edge of the theoretical uncertainty band. One can thus 
conclude that no serious disagreement is present, but that the situation certainly 
deserves further investigations, also in the light of data on forward bottom pro- 
duction by the DO Collaboration, which appear to be about a factor of four higher 
than the theoretical prediction. 



II RESUMMATION OF THRESHOLD EFFECTS 

As we mentioned in the Introduction, potentially large logarithms appear in the 
NLO fixed-order calculation and make it less reliable in some regimes. We deal 
now with log(l — 4m 2 /s) terms, which appear in the partonic cross section in the 
form 



a(s) = a (s) 



l + CYUog 2 (l-— ) + 



(2) 



These logarithms become large when the cms energy the colliding partons have is 
close to the invariant mass of the two produced heavy quarks, that is to say, when 
the quarks are produced close to the threshold. Resummation of such terms can 
turn out to be phenomenologically important when the quark is very heavy with 
respect to the hadronic collision energy: top production at the Tevatron, with the 
top weighing 175 GeV and the machine delivering 1800 GeV, could be such a case. 

In order to better investigate this issue three groups have in recent years at- 
tempted such a resummation. We are not going to go into any details here, to be 
found in the respective papers, but simply summarize their conclusions and point 
to similarities and differences of approaches and results. 



Laenen, Smith and van Neerven [11] performed the resummation in x-space, 
by directly exponentiating the leading log (LL) a s log 2 (l — x) term, where x = 
4m 2 /s. Their formulae necessitate the introduction of an artificial infrared cutoff 
which serves a double purpose: it avoids hitting the Landau pole when the now 
z-dependent argument of the running coupling a s tends to zero as z — > 1, and it 
regulates an otherwise divergent integral of the form Jq dzexp[\a\ log 2 (l — z)\. They 
predict a 10% increase over the NLO fixed order calculation for top production at 
the Tevatron, but with a large uncertainty due to the dependence on the unphysical 
cutoff. 

Berger and Contopanagos [12] instead perform the resummation in Mellin mo- 
ments space, and then invert back to x-space for producing phenomenological pre- 
dictions. When performing this inversion they apply a so-called principal value 
prescription to deform the integration contour and avoid hitting the Landau pole. 
They then discard all the non-leading terms generated by the inverse transform, 
arguing they are not universal, and only retain the LL ones. The further need for an 
infrared cutoff to avoid the Jq dzexp[\a\ log 2 (l — z)\ divergence is met by choosing it 
in such a way that all integrations are confined to a perturbative domain, i.e. to re- 
gions where the discarded non-leading terms are really subdominant. Through this 
procedure they also find an increase of about 10%, but claim a small uncertainty 
due to the motivated choice of the cutoff. 

A third evaluation of soft-gluon resummation effects has been performed by 
Catani, Mangano, Nason and Trentadue [13]. They also perform the resummation 
in Mellin space, and avoid the Landau pole by what they call minimal prescription, 
i.e. a choice of the integration contour which leaves the non-perturbative branch 
cut to the right. After doing so, they perform the z integration by retaining both 
the leading and the next-to-leading (NLL) terms. Their argument is that the NLL 
contributions are generated by momentum conservation and, while formally sub- 
leading, discarding them leads to factorially divergent integrals and hence to the 
need for a cutoff. The phenomenological outcome of their investigation, at vari- 
ance with previous findings, is that top hadroproduction at the Tevatron energy is 
predicted to increase by only 1% with respect to the fixed order result. They also 
claim little uncertainty on this result due to the absence of an explicit cutoff. 

Whatever the approach one considers correct, we are however still far from ex- 
perimentally probing such fine details (see experimental errors in fig. 1). 

Ill LARGE TRANSVERSE MOMENTUM 

RESUMMATION 

Another kind of large logarithmic terms appearing in the NLO calculation and 
eventually spoiling its convergence are \og(p^/m 2 ) terms. They also need therefore 
to be resummed to all orders to allow for a sensible phenomenological prediction. 
Such a resummation has been performed along the following lines [14] . 
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FIGURE 4. Comparison of bottom pr distributions in the fixed order (NDE) and resummed 
approaches (left plot, from ref. [14]), and effect of a harder b — > B fragmentation function on 
bottom production at the Tevatron [18] (right plot). 

One observes that in the large-pr limit {jpr ^> m) the only important mass terms 
are those appearing in the logs, all the others being power suppressed. This means 
that an alternative description of heavy quark production can be achieved by con- 
sidering massless quarks and providing at the same time perturbative distribution 
and fragmentation functions also for the heavy quark, describing the logarithmic 
mass dependence. The factorization formula becomes 

d<r(p T ) = H j F i/uMi [™\) F j/H*(n, [m])dfr(ij -> k;p T , /x)D?(/x,m), (3) 

ijk 

with parton indices i,j and k also running on Q, taken massless in a, now an MS 
subtracted cross section for light parton production. The dependence on m of the 
parton distribution functions Fi/n, shown between square brackets in eq. (3), is 
only there when i or j happens to be the heavy quark Q. 

The key point is that the large mass of the heavy quark allows for the evaluation in 
perturbative QCD (pQCD) of its distribution and fragmentation functions. Initial 
state conditions for Fq/ h {^,q = m) [15] and D^(/i ~ m) [16] can be calculated 
in pQCD at NLO level in the MS scheme. These initial state conditions can then 
be evolved with the Altarelli-Parisi equations up to the large scale set by fx ~ pr- 
This evolution will resum to all orders the large logarithms previously mentioned. 

It is important to mention that due to the neglecting of power suppressed mass 
terms this approach becomes unreliable when — m. In this region only a case 
by case comparison with the full NLO massive calculation - here reliable and to be 
taken as a benchmark - can tell how accurate the resummed result is. 

Phenomenological analyses [14] show that the effect of the resummation becomes 
sizeable only at very large px, say greater than 50 GeV for bottom hadropro duct ion 



at the Tevatron, resulting in a milder factorization/renormalization scale depen- 
dence of the result and in a slightly softer p T spectrum. 

Non-perturbative effects describing, say, the b — > B meson transition can also 
be included within this formalism. The b — > B fragmentation function can be 
fitted to LEP e + e~ data and then used for predicting B cross sections at the 
Tevatron. Recent analyses [17] show that this FF, when used in connection with a 
NLO evaluation of bottom production like the MNR one, should probably be taken 
harder than commonly done in the past. This choice of a harder non-perturbative 
FF would increase the Tevatron cross sections from twenty to fifty per cent [18] 
(see fig. 4), helping to reconcile theory and experimental data. 

IV HEAVY QUARKONIUM PRODUCTION 

The production of heavy quarkonia has been subjected to intense study in the 
last two or three years, with tens of papers having being produced on the problem 
of J/ijj, ip', x an d T production in e + e~ , jp, pp, pN, tiN collisions. 

The reason for such a surge in interest was the appearance of a theoretical frame- 
work, the so called Factorization Approach (FA) by Bodwin, Braaten and Lepage 
[19], which seems able to solve the theoretical problems that quarkonium produc- 
tion models faced in the past, and also to reconcile theoretical predictions with 
experimental data, previously in disagreement up to factors of fifty in some in- 
stances. 

In this talk I shall not review the Factorization Approach in detail, leaving this 
theoretical introduction to other sources (see for instance [20]). I shall just recall 
how the FA writes the quarkonium state H production cross section in the following 
form: 

a{ij -> QQ _> H) = 5>(ij - QQ[n])(0 H (n)). (4) 

n 

According to this equation, the cross section for producing the observable quarko- 
nium state H is factorized into two parts. In the short distance part a QQ pair of 
heavy quarks is produced in the spin/colour /angular momentum state 2S+1 L^ = n 
by the scattering of the two light partons % and j. Subsequently this pair hadronizes 
into the quarkonium H; (0 H (n)) is formally a non-relativistic QCD (NRQCD) ma- 
trix element describing this non-perturbative transition. 

An important feature of this equation is that also QQ pairs in a colour octet 
state are allowed to contribute to the production of a colour singlet quarkonium H: 
their colour is neutralized via a non-perturbative emission of soft gluons. While 
the corresponding matrix elements are suppressed by the need of such an emission, 
the short distance coefficients can on the other hand be large, perhaps overcom- 
pensating the suppression of the non-perturbative term. This explains why colour 
octet contributions can play a very important role in predicting the total size of 
quarkonium production cross sections. 



The limited space available doesn't allow any detail about the phenomenolog- 
ical studies which have been performed. Ref. [21] contains a small up-to-date 
review, with references to the original papers. I shall only mention here that the 
introduction of the colour octet channels allows for what looks like a "reasonable" 
description of the data. Colour octet matrix elements have to be fitted to the data 
themselves, and the uncertainties on these fits are certainly not smaller than a 
factor of two, due to the many systematics entering their determination: parton 
distribution functions, charm quark mass, a s value, higher order QCD corrections, 
etc. However, the values one gets appear of the correct order of magnitude if com- 
pared via NRQCD scaling rules [22] to the colour singlet ones, thereby supporting 
the underlying picture. 

So far all fits have been performed using leading order cross sections. Very 
recently a next-to-leading order calculation for quarkonium total cross sections, 
both via singlet and octet channels, has been completed [23]. It will therefore 
be possible to reduce at least some of the uncertainties previously mentioned and 
hence obtain more reliable fits. 
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